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Abstract — Shadowing is believed to degrade the quality of 
service (QoS) in wireless cellular networks. Assuming log-normal 
shadowing, and studying mobile's path-loss with respect to the 
serving base station (BS) and the corresponding interference 
factor (the ratio of the sum of the path-gains form interfering 
BS's to the path-gain from the serving BS), which are two key 
ingredients of the analysis and design of the cellular networks, 
we discovered a more subtle reality. We observe, as commonly 
expected, that a strong variance of the shadowing increases 
the mean path-loss with respect to the serving BS, which in 
consequence, may compromise QoS. However, in some cases, 
an increase of the variance of the shadowing can significantly 
reduce the mean interference factor and, in consequence, improve 
some QoS metrics in interference limited systems, provided the 
handover policy selects the BS with the smallest path loss as the 
serving one. We exemplify this phenomenon, similar to stochastic 
resonance and related to the "single big jump principle" of the 
heavy-tailed log-nomal distribution, studying the blocking prob- 
ability in regular, hexagonal networks in a semi-analytic manner, 
using a spatial version of the Erlang's loss formula combined with 
Kaufman-Roberts algorithm. More detailed probabilistic analysis 
explains that increasing variance of the log-normal shadowing 
amplifies the ratio between the strongest signal and all other 
signals thus reducing the interference. The above observations 
might shed new light, in particular on the design of indoor 
communication scenarios. 

Index Terms — Wireless cellular networks, blocking probabil- 
ity, path-loss, shadowing, indoors, interference factor, stochastic 
resonance, geometry, honeycomb, Poisson. 



I. Introduction 

Modeling of the attenuation of an electromagnetic wave as 
it propagates in space is a major component in the analysis 
and design of wireless systems. This phenomenon, also called 
propagation loss, is caused by the decay of the signal power 
with the distance from the emitter (existing even in the 
free space propagation models) and due to various obstacles 
between emitters and receivers (trees, buildings, hills, etc.) 
present in real network profiles. Complexity and haphazard 
character of actual network profiles makes pertinent the sta- 
tistical modeling of the propagation loss. In this approach, 
the propagation loss between an emitter and a receiver, called 
path-loss, is typically modeled by the product of the distance- 
loss function — a deterministic function of the distance 
between the two antennas, which represents average path-loss 
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on the given distance in the network, and a random variable, 
called shadowing, that takes into account in a statistical 
manner the deviation from this average, observed for each 
particular pair of emitter and receiver We call this model path- 
loss with shadowing. The distance-loss function is commonly 
assumed to be some power of the distance, with the expo- 
nent called path-loss exponent. The random variable of the 
shadowing is often assumed to have log-normal distribution, 
normalized to have mean one and parametrized by its variance 
or standard deviation. 

Various QoS metrics in cellular networks, as blocking prob- 
ability for constant bit-rate (CBR) connections and spectral 
efficiency for variable bit-rate (VBR) connections, depend on 
the strength (i.e., variance) of the shadowing. It is commonly 
believed that an increase of the variance of shadowing pe- 
nalizes the network performance. The results presented in 
this paper shed some new light on this problem. Namely, 
studying the blocking probability (defined as the fraction of the 
CBR connections that cannot be established due to insufficient 
transmission resources, in the long run of the system) we 
have discovered that it is not always increasing with the 
variance of the shadowing. For example, in our model of 
the OFDMA hexagonal network consisting of 36 BS, with 
cell radius 0.525km and the path-loss exponent equal to 2.5, 
the blocking probability evaluated at the presence of the log- 
normal shadowing with the standard deviation of 25dB is four 
times smaller than in the scenario with no-shadowing. Even 
if this spectacular example regards a very strong shadowing, 
we obtain a smaller, but still very significant, decrease of 
the blocking probability for the shadowing with the standard 
deviation from 7 to 15dB, which might be appropriate for the 
indoor scenario (user-indoors, BS-outdoors); cf [ ]. 

In all cases, a very strong shadowing ultimately makes the 
blocking probability tend to 1 and this dependence indeed 
becomes (as expected) monotone for higher path-loss exponent 
(larger than 4 in the considered examples). 

To explain the above, somewhat surprising, observations and 
extend them to other QoS metrics we study the impact of the 
shadowing and the path-loss exponent on the following two 
key characteristics of any given mobile in the network: 

• its path-loss to the serving BS, which is the one received 
with the strongest signal (and not necessarily the closest 
one), 

• the so called mobile's interference factor, defined as the 
ratio of the sum of the path-gains form interfering BS to 
the path-gain from the serving BS. 

These are two key ingredients in the analysis of wireless 
cellular networks and thus their mean values can be con- 



sidered as some QoS "pre-metrics". In particular, they are 
explicitly present in the call blocking condition — the one 
used to control the admission of streaming users, and hence 
intrinsically related to the blocking probability. They are, even 
more straightforwardly, related to the spectral efficiency of 
the data networks. While being key ingredients in the study 
of wireless systems, the (mean) QoS pre-metric are also more 
easy to analyse. In particular they do not depend on a particular 
assumption regarding the spatial correlation of shadowing 
(which is not the case, e.g. for the mean blocking probability). 

We have studied the mean values of the above two basic 
QoS ingredients, with the averaging taken over all possible 
locations of users in the network and over the distribution of 
the shadowing. 

Our main findings are as follows: 

• The mean path-loss (with respect to the serving BS) is 
always increasing in the variance of the shadowing. The 
ultimate degradation of the QoS for large shadowing vari- 
ance is due to this increasing path-loss. (When possible, 
this may be however remedied by increasing the power 
of the emitted signals). 

• The mean interference factor is not monotonic in the 
variance of the shadowing. It first increases and then 
decreases (asymptotically to zero!), when the shadowing 
variance goes to infinity. This asymptotic behaviour can 
be heuristically explained by the single big jump principle 
of heavy-tailed distributions: the sum of the (log-normal) 
path-gains form all antennas is dominated by a big value 
of the path-gain from a single antenna. When this antenna 
is the serving one, then this big path-gain does not 
count in the interference, which becomes negligible in 
proportion to it. 

• The above two facts lead to the phenomenon that we 
may call a stochastic resonance for QoS in path-loss-and- 
interference limited systems: when QoS is not yet com- 
promised by path-loss conditions, a moderate increase 
of the shadowing variance may make it profit from the 
reduction of the interference. 

We confirm the above findings by a mathematical analysis 
of the respective stochastic models. We also compare in this 
matter the performance of the perfect (hexagonal) and irregular 
(Poisson) networks and find that both architectures exhibit 
very similar QoS "pre-metrics" for the standard deviation 
of the shadowing larger than 20dB. Moreover, we prove 
an interesting invariance of the QoS metrics of the infinite 
Poisson cellular networks with respect to the distribution of 
the shadowing. As a consequence we also obtain fully explicit, 
analytical results for the mean path-loss and interference 
factors in the case of the infinite Poisson network. 

The remaining part of this paper is organized as follows. 
In the next section we briefly present related works. In 
Section III we describe our models. The main numerical 
results are presented in Section IV. Next, in Section V we 
present mathematical analysis of the models, which supports 
and completes our numerical findings. Finally, in Section VI 
we provide some concluding remarks. 



II. Related works 

The propagation loss model considered in this paper is com- 
monly accepted in the literature; see e.g. [ ] where log-normal 
shadowing of mean 1 is considered. A possible extension of 
this model consists in assuming shadowing distribution (say, 
its variance) that depends on the distance, cf. [ ]. 

The impact of the shadowing on the distribution of the 
interference factor is studied numerically in [ ] and analyt- 
ically in [.,]. However, the above two articles do not take into 
account the modification of the network geometry induced by 
the shadowing, i.e., assume that mobiles are served by their 
geographically closest BS. This is not a realistic assumption 
and, as we will show in this paper, leads to misleading 
conclusions that the shadowing dramatically increases the 
mean interference factor. 

The paper [7] focuses on the interference factor averaged 
over a given cell, and in particular the effect of shadow- 
ing on this average. It is shown there that the cell shape 
modification induced by the shadowing affects significantly 
the mean interference factor. More precisely, that this mean 
decreases substantially if mobiles are served by the BS offering 
the smallest path-loss. We adopt this assumption throughout 
the present paper in the context of regular (hexagonal) and 
irregular (Poisson) geometry of BS, as proposed in [ ]. 

Some papers (see e.g. [9, 10]) propose more explicit approx- 
imations of the interference factor and its moments (mean and 
variance) assuming only deterministic propagation loss models 
(without random shadowing). [ ] studies the distribution of 
the interference factor in such a case. 

In [ ] the authors partially confirm, by a different ap- 
proach, our early observation from [1], that the average SIR 
(which is the inverse of the interference factor) might increase 
with the shadowing variance when the best server policy is 
chosen. 

The interference factor was recognized very early as a key 
element in the performance evaluation of cellular networks; 
cf. [13, 14]. Fundamental to our approach to the evaluation 
of the blocking probability are papers [15, 16]. They show 
how the power allocation problem without power limitations 
can be reduced to an algebraic system of linear inequal- 
ities. Moreover, they recognize that the spectral radius of 
the (non-negative) matrix corresponding to this system not 
greater than 1 is the necessary and sufficient condition of 
the feasibility of power allocation without power limitations. 
This approach lead to the development of a comprehensive 
framework of the evaluation of the blocking probability in 
CDMA, HSDPA and OFDMA, via a spatial version of the 
famous Erlang's formula in [8, 17-19], QoS in data networks 
are studied using this approach in [20]. 

Finally, recalling that the mean QoS pre-metrics studied 
in this paper do not depend on the spatial correlation of the 
shadowing, we remind [2 1 , 22] as bringing models that can be 
used when studying the spatial distribution of the QoS metrics. 

III. Model description 

A. Location of base stations 

In this paper we will consider two particular models for 
the location of BS, hexagonal and Poisson one. The former 



is commonly considered as an "ideal" model for the cellular 
networks, while the latter one can be seen as an extremal case 
of very irregular network. 

1} Infinite Models: 

• Hexagonal network. Consider BS located on a regular 
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Fig. 1. Hexagonal pattern of 4 X 4 BS on rectangular torus T4. Identified 
points are denoted by the same digits or characters. 

hexagonal grid on M^ with the distance A between two 
adjacent vertices of this grid '; cf. Figure 1. Note that 
the surface area of a given cell (hexagon, i.e. subset of 
the plane whose points are closer to a given point of 
the grid than to any other) of this model is equal to 
\/3A^/2. Thus the intensity of the BS in this model 
is equal to A = 2/(^/3A2) BS/km^. In what follows it 
will be customary to consider a stationary version $/f 
of this grid, which can be obtained by randomly shifting 
it through a vector uniformly distributed in one given 
hexagon (cf. [23, Example 4.2.5]). In this model a given 
location, say the origin of the plane, corresponds to an 
"arbitrary" location of a mobile, "randomly chosen" in 
the network. 
• Poisson network. Assume that BS are located at the 
points of a stationary, homogeneous Poisson point process 
(p.p.) $p of intensity A BS/km^ on the plane M^. 
When comparing performance of Poisson and hexagonal 
model we will always take them with the same intensity 
A = 2/(^A2). 

Considering infinite models is often a convenient way of 
studying phenomena arising in very large networks. A particu- 
lar property of these models is lack of (geographic) boundary 
effects, which in real, large but finite, networks, have often 
a negligible impact on performance characteristics measured 
in the "middle" of the network. However, as we will see in 
this paper, sometimes mathematical assumption of an infinite 
network may create some artifacts, which are not observed in 
more realistic, large but finite, networks. 

2) Bounded Models: In order to have finite network mod- 
els, and still neglect the boundary effects (which might be rea- 
sonable for large networks) one often considers toroidal model. 

' The set of vertexes of this grid can be described on the complex plane by 



Recall that, roughly speaking, rectangular torus is a rectangle 
whose opposite sites are "identified". For iV = 2,4,6,..., 
we will denote by Tj^ the rectangle [-A^A/2,iVA/2) x 
[-iV\/3A/4, iV\/3A/4) with toroidal metric. Restricting ^n 
to Tjv, i-C- taking ^j^ ~ $/f n Tat one obtains the model 
whose distribution is invariant with respect to translations on 
the torus. Thus we obtain a hexagonal network model that 
consists of N"^ cells (cf. Figure 1) and which does not exhibit 
any border effects. Similarly we will consider the restriction 
$p of the Poisson p.p. $p to Tjv- 

B. Path-loss model with shadowing 

For a given BS X e $ ($ = $p or (^h) and a given 
location y e M^ on the plane we denote by Lx{y) the (time- 
average, i.e., averaged out over the fading) path-loss between 
BS X and location y. In what follows we will always assume 
that 
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where £(•), called distance-loss, is a non-decreasing, de- 
terministic function of the distance between an emitter and 
a receiver, and S'x(-) is a random shadowing field related 
to the BS X. In what follows we call Lx{y) path-loss 
with shadowing (or path-loss for short) between X and y. 
Moreover, we will always assume that given locations of BS 
{Xi S $} their shadowing fields {Sxi{-)} are independent 
non-negative stochastic processes, each being indexed by 
locations y S M^. More formally speaking, the locations of 
BS X and their respective shadowing fields S'x(-) form an 
independently marked version $ = {(X, S'x(-))}xe* of the 
point process $. 

Regarding the distribution of the marks (shadowing fields) 
of this process, they are assumed to have the same marginal 
distributions, i.e., given X, Sx{y) has the same distribution 
for all y e M^, of normalized mean E[S'x(2/)] = 1> with the 
following two cases being of particular interest 

• Sx{y) = 1, which corresponds to a case with negligible 
shadowing (we will say also "no shadowing"), 

• for all y, Sx{y) is log-nornal random variable with 
mean 1. Recall that such a mean-1 log-normal variable 
S can be expressed as 5 = 6^+°"^ where N is standard 
Gaussian random variable (with mean and variance 1) 
with p = — (7^/2 and some constant a. Indeed, in this 
case E[S] = gAi+^Vs ^ i Note that if the shadowing 
is log-normal random variable then the path-loss (at a 
given distance) expressed in dB is Gaussian random 
variable. Furthermore, in this context it is common to 
parametrize the log-normal shadowing by the standard 
deviation (SD) of S expressed in dB, i.e., the SD of 
10 logj^Q S. We will denote it by v. With respect to 
the previous parametrization we have v = crlO/loglO. 
Throughout the paper we will call v the logarithmic 
standard deviation (log-SD) of the shadowing. 

If not otherwise specified, we do not make any particular 
assumption on the correlation of the shadowing field Sxiy) 
for given X and different locations y. 



Throughout the paper we will implicitly assume also that 
mean path-gain is finite, i.e., E[l/S'] < oo. Note that this 
condition is satisfied for log-normal variable, indeed, in our 
case of mean-1 variable E[l/S'] = e'"' = e''^^°s^ w/wa^ 

For the deterministic distance-loss function L() the fol- 
lowing particular model is often used and will be our default 
assumption in this paper: 

L{r)^{Krf (2) 

where K > and /3 > 2 are some constants. 



C. Handover policy and path-loss factor 

In what follows we will assume that each given location 
y e M^ is served by the BS X* e $ with respect to which 
it has the weakest path-loss Lx*{y) (so, in other words, the 
strongest received signal, given all BS emit with the same 
power), i.e, such that 



Lx'iy)<Lx{y) 



for all X e $ . 



(3) 



with any tie-breaking rule. Note that in the case of negligible 
shadowing {Sx{y) = 1) and strictly increasing function L(-) 
the above policy corresponds to the geographically closest BS. 
Note also that for our infinite network models with random 
shadowing, one has to prove that the minimum of the path-loss 
is achieved for some BS, i.e., that X* is well defined. 

Note that Lx' (y) is the path-loss experienced by a user 
located at y with respect to its serving BS. Obviously it 
determines the QoS of this user (we will be more specific 
on this in Section III-E). In this context we will call it path- 
loss factor- of user y and denote by l{y) = Lx*{y)- Note 
that it depends on the location y but also on the path-loss 
conditions of this location with respect to all BS in the network 
l{y) = l{y, $). Path-loss factor l{y) is typically not enough to 
determine the QoS of a given user. 



D. Interference factor 

For a given location y E 
factor f{y) as 



we define the interference 
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(4) 
provided X* is well defined. 

Study of the path-loss and interference factors, which are 
relatively simple objects, can give an important insight into 
more involved QoS metrics, such as blocking probability in 
streaming traffic and mean throughput in data traffic. In what 
follows we recall how l{y) and f{y) appear naturally in the 
evaluation of the blocking probabihties. 

^not to be confused with the path-loss exponent /3 



E. Blocking probability; a space-time scenario 

In this section we briefly describe the relation between the 
path-loss and interference factors and the blocking probability. 
This relation, whose very essence can be explained by the 
famous Erlang's loss formula, was observed in the current 
geometric context (however without shadowing) in [ ]. 

In order to evaluate the blocking probability it is necessary 
to specify the dynamics of call arrivals and their durations, 
as well as to identify the set of feasible configurations of 
users (which can be served simultaneously at their requested 
bit-rates). To this regard, consider a given realization of the 
netowrk with shadowing $, and a spatio-temporal Poisson 
arrival process of calls which require from the network some 
predefined transmission rates for some exponential transmis- 
sion times. These rates can be maintained at the price of 
blocking of some call arrivals when a network congestion 
occurs. The fractions b = 6(<i>) of blocked arrivals in the long 
run of the system is called the blocking probability. By the 
famous Erlang 's loss formula, it is equal to the conditional 
probability that in the stationary configuration of the (non- 
blocked) arrival process the system cannot admit a new user, 
given all users in the current configuration can be served. 
Moreover, if the decision whether to block a given call (or 
admit it) is based on the verification of some feasibility condi- 
tion that has the so called multi-Erlang form, then the Erlang's 
loss formula can be relatively easily evaluated, e.g. discretizing 
the values of the SINR and using Kaufman-Roberts algorithm. 
A canonical form of the multi-Erlang feasibility condition 
involves verification by each BS X of the following condition 



v-x:=x 



{Ky)Jiy) 



< 1, 



(5) 



where the summation is taken over all users (including a 
new arrival) to be served by the BS X and (/?(•,•) is some 
function of the path-loss and interference factors of user 
y. This condition guarantees sufficient wireless resources to 
maintain the predefined transmissions rates for all served 
mobiles. Specific form of the function (p{-,-) needs to be 
developed for each particular cellular technology (taking into 
account the performance of the coding schemes, type of the 
multiplexing, etc.). Below we show two examples borrowed 
from our previous studies. They give some insight into how 
the feasibility condition (5) depends on the user transmissions 
rates, it is supposed to guarantee. 



• For the down-link in the OFDMA network 



filJ) 



W^j{{l-e)/{im/P) + a + f) 



(6) 



where P is the maximal BS power, e is the fraction of 
this maximal power used in common (pilot) channels, 
a is the intra-cell orthogonality factor (usually assumed 
to be in OFDMA), N external noise power, W is the 
system bandwidth, r is the required bit-rate r of user and 
V' is the link performance function (V'(^) is the bit-rate 



per Hz available when SINR is equal to ^; ^y, cf. [19]. 
For the down-link in CDMA network 



^(IJ) 



/Nl 



l + a^l-e\ p 



^o./) 



(7) 



where ^ — ip~^{r/W) is the SINR threshold correspond- 
ing to the required bit-rate r of user and the remaining 
notation notation as above; (cf. [18]). 

In what follows we will denote by E[&] = E[6($)] the 
blocking probability averaged over possible scenarios re- 
garding locations of BS and their shadowing conditions. It 
can be evaluated by the simulation of several realizations 
of the network with shadowing $, evaluation of 6($) by 
the Kaufman-Roberts algorithm as described above, and then 
taking the empirical average over the realizations of $. How- 
ever, in practice one realization of $ is enough, provided the 
shadowing fields Sx (y) do not exhibit high spatial correlation 
across y ( recall that we have assumed them to be independent 
across X); cf. Footnote 4. Indeed, we have noticed in our 
experiments, that for large enough networks (in the case of the 
hexagonal network Tg is enough!) with spatially uncorrected 
shadowing, the value of 6($) is almost invariant with respect 
to $ and hence very close to E[6(<&)]. This is due to spatial 
ergodic properties of the process $. 



F. Our methodology in the study of the network QoS 

In section IV-A we will show some numerical examples, 
which show the typical dependence of the blocking probability 
E[fo] on the parameters of the path-loss model. These exam- 
ples, restricted to OFDMA, are not supposed to be exhaustive. 
The goal is to show the typical tendencies. 

In order to explain these tendencies, in Sections IV-B 
and IV-C we will study more thoroughly the mean values of 
the interference and path-loss factor E[/(2/)] — E[/(y, $)], 
E[/(y)] = E[;(y, $)] (where the expectation E[. . .] coiTe- 
sponds to the distribution of $, i.e., this of the shadowing 
field and of the random location of the user). By the translation 
invariance of the distribution of our infinite and toroidal mod- 
els, these expectations (corresponding the spatial averaging) 
do not depend on the user location and thus, for these models. 



e.g., assuming additive wliite Gaussian noise (AWGN) cliannel and tlie 
link performance closed to the optimal one, ip is given by the famous 
Shannon's formula tp{^) = log2(l + ^). Taking V(C) = 'ilog2(l + 5) 
with some constant a < 1 permits to account for a degradation of the 
link performance in practical systems compared to the ideal AWGN case; 
cf. [24]. Further extensions consider the Single-Input-Single-Output (SISO) 
AWGN channel with fading, for which the known formula for the ergodic 
capacity is </>(§) = E[log2(l + C|i^P)], where the expectation is with 
respect to the distribution of the channel fading F, and the Multiple- 
Input-Multiple-Output (MIMO) AWGN channel, whose ergodic capacity is 
tp{£,) = E[log2 det(/ + gFF-^)], where F is the vector of channel fading; 
cf. [25]. 



nKy)] = nm] and E[/(y)] = E[/(0)] \ 

Our methodological conjecture is as follows. We believe that 
the mean path-loss and interference factors E[/(0)], E[/(0)] 
can be considered as primitive (basic) metrics of the QoS 
and their behavior can (at least qualitatively) explain the 
main tendencies observed for more involved QoS metrics. This 
methodological conjecture is motivated by the observation that 
the function (^ in the feasibility condition (5) is an increasing 
function of some linear combination of l{y) and f{y) (at 
least for the examples of CDMA and OFDMA given above). 
Indeed, we will show that the study of E[?(0)] and E[/(0)] can 
explain the aforementioned non-monotonicity of the blocking 
probability E[6] with respect to the standard deviation of the 
shadowing. 

IV. Numerical results 

Following the methodology described in Section III-F, in 
this section we will first study the blocking probability and 
then the mean path-loss and interference factors for the hexag- 
onal and Poisson network models with log-normal shadowing. 

A. Blocking probability 

In this section we consider only the hexagonal network on 
the torus Tg. We evaluate the blocking probability E[&] in 
OFDMA network using the Kaufman-Roberts algorithm, as 
described in Section III-E. 

We assume the following parameter values for OFDMA: 
System bandwidth W — 5MHz. BS are equipped with omni- 
directional antennas having a gain 9dBi and transmit with the 
maximal power 43dBm; thus P = 43 + 9 = 52dBm when we 
account for the isotropic antenna gain. The common channel 
power is the fraction e = 0.12 of P. The ambient noise power 
is assumed N ~ — lOSdBm. 

We assume perfect intra-cell orthogonality, i.e., a = 0. 
Moreover we are not considering any opportunistic scheduling 
over fading. This allows us to characterize the link per- 
formance (averaged over fading) via the Shanon's formula 
r/W = -0 [0 — log2 (1 + 0' without specifying into how 
many sub-carriers the bandwidth W = 5MHz is split. 

We assume a traffic demand of 46.2 Erlang per km^ 
consisting of streaming calls at the bit-rate r = 180Kbits/s 
(typical for videoconferencing) that is served by the hexagonal 
network consisting of 36 BS (on the tours Tg) with the distance 
between adjacent BS A = 1km. 

The (deterministic) distance-loss function is L{x) = (Kx) 
with K — 8667km~^ (which follows from Cost-Hata 

"* Often the mathematical expectation E[/(0, #)] (and similarly 
for E[/(0, $)]) corresponds to the empirical mean value 
lim„— >oo l/iX] /(S/ii '^) °f 'he interference factor measured at many 
locations "uniformly" sampled in one given realization of the network and 
shadowing. A precise statement and rigorous proof of such an ergodic result 
is beyond the scope of this paper. We remark only that for the hexagon 
network on the torus, this result follows simply form the Law of Large 
Numbers, when j/i ai'e independently and uniformly distributed and provided 
the shadowing variables SxiVi) are independent across different values of 
yi. Indeed, in this case f{yi,^) and l{yi,^) are independent, identically 
distributed (across i) random variables. However, recall that the latter 
assumption, corresponding to spatially uncorrelated shadowing, is not our 
default assumption, since it is not needed for other results regarding E[J(0)] 
and E[/(0)]. 




Logarithmic standard deviation of the shadowing 

Fig. 2. Blocking probability in OFDMA hexagonal network on the torus Tg 
with log-normal shadowing with log-SD v and path-loss exponent /3, evaluated 
using the Kaufman-Roberts algorithm for the traffic 46.2 Erlang per km^ . 



model [26] for urban areas, assuming frequency 1795Mhz, 
BS antenna height 50m, mobile antenna height 1.5m, for /3 = 
3.38). Moreover, we assume that the values of the shadowing 
Sx {y) for given X and different locations y are independent. 
Figure 2 shows the dependence of the blocking probability 
E[fo] on the path-loss exponent /? and the logarithmic standard 
deviation v of the log-normal shadowing evaluated in our 
OFDMA network model. 

Remark 4.1: • For negligible shadowing (logarithmic 
standard deviation close to 0) the blocking probability 
E[6] first decreases in the path-loss exponent ji (on our 
figures for /3 from 2.5 to 4) and then increases in /?. 
• The blocking probability is not always increasing in 
the standard deviation of the shadowing. Indeed, on our 
figures with /3 > 4 it is monotone increasing. However, 
for P < 3.5 the blocking probability E[6] first increases, 
then decreases, and ultimately increases to 1. 
Note the decrease of the blocking probability in the standard 
deviation of the shadowing can be quite significant even 
between 7 and 15 dB, depending on the path-loss exponent. 

The lack of monotonicity observed in Remark 4.1 is not 
specific for our choice of the traffic of 46.2 Erlang par km^ 
as can be remarked on Figure 3, where we have assumed 
two different smaller values of the traffic. Moreover, we have 
observed very similar patterns, not presented here due to 
space constraints, in our model of CDMA. Furthemore, we 
have confirmed these results by the the crude Monte-Carlo 
simulations of the network with the arrivals and departures of 
users (implemented in MATLAB). 

In the next section we will explain this behavior and argue 
that it my be expected for other QoS metrics which depend on 
some combination of the path-loss factor and the interference 
factor. 

B. Analysis of the interference factor 

Now, we will study the impact of the shadowing and also 
the geometry and size of the network on the interference factor 
E[/(0)] that is a key to the understanding of the strange non- 
monotonicity of the blocking probability shown above. Recall 
that, contrarily to the blocking probability, the expectation 
E[/(0)] (as well as E[^0)]) does not depend on any particular 



correlation of the values of the shadowing Sx (y) for given X 
and different locations y. 

We begin with an important observation made directly from 
our model. 

Remark 4.2: By the homothetic invariance of our hexago- 
nal and Poisson models on the torus, or in the infinite models, 
with the distance-loss function (2), the mean interference 
factor does not depend on the intensity A of BS but only on 
the size N of the network. 

Figures 4(a) and 4(b) show the impact of the path-loss 
exponent, shadowing and the size of the network in the case of 
the hexagonal and Poisson network architecture, respectively. 
Here are our main observations. 

Remark 4.3: 1) Observe on Figure 4(a) for hexagonal 
network of a given size N^ BS, with A^ = 6, 10, 30, 
and a given path-loss exponent /3 = 3,4, 5, that the 
mean interference factor 'Ei[f (0)] first increases and then 
decreases to when the value v of logarithmic standard 
deviation (log-SD) of the shadowing increases. 

2) For the Poisson network (see Figure 4(b)) E[/(0)] 
decreases in log-SD starting already from very small 
values of v. 

3) The actual size of the network consisting of N^ BS, when 
N > 100, has negligible impact on E[/(0)] when j3 — A 
and u < 10 or /3 = 5 and u < 15 both in hexagonal 
and Poisson case (in this latter case N"^ is the expected 
number of BS). In this regime the value of E[/(0)] 
corresponds to this in the respective infinite model. In 
particular, for Poisson network it is equal to 2/(/3 — 2) 
and does not depend on log-SD v (cf. Proposition 5.6 
below). 

4) When (3 — 4: and v > 10 or /3 = 5 and w > 15 the mean 
interference factor E[/(0)] non-negligibly increases with 
the network size. 

5) Comparing Figures 4(a) and 4(b) for w > 20 we 
observe that /or large log-SD of the shadowing the mean 
interference factor evaluated for the Poisson network is 
almost exactly the same as for the hexagonal network 
of the same size. 

Remark 4.4: The seminal paper [ ] considers only the 
hexagonal network architecture, however, the beneficial impact 
of the shadowing is not observed there. The reason is that the 
model considered in [ ] assumes that the smallest-path-loss 
BS (the serving one) is selected among the Nc closest BS. 
In particular, Nc = 1 ignores the shadowing in the handover 
policy as it corresponds to the situation where the serving BS is 
always the closest one. On the other hand the model considered 
in our paper corresponds to Nc equal to the total number of 
BS in the network. In consequence, for a higher path-loss 
exponent (say /3 = 4) and small and moderate log-SD of the 
shadowing (0 < w < 12) our numerical results are close to 
those of [ ] with Nc = 4; cf. our Figure 4(a) and the last 
column in Table 1 in [;]. The fact that the average interference 
factor decreases in some cases with log-SD of the shadowing 
has not been observed in [7] due to the set of parameters 
considered there. Indeed, for a smaller path-loss exponent, 
(3 — 3, our Figure 4(a) shows the mean interference factor 
decreasing in v starting from u w 8. This range of parameters 
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Fig. 3. Blocking probability for OFDMA network as on Figure 2 with traffic 34.6 Erlang per km^ and 23.1 Erlang per km^ 
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Fig. 4. Mean interference factor in hexagonal and Poisson network on the torus Tjv with log-normal shadowing with log-SD v and path-loss exponent /3. 
Note that E[/(0)] increases with the size of the network. The straight lines correspond to the infinite (on R^) Poisson model; cf. Proposition 5.6. 



is also considered in [ , Table 2] however, with the Nc = 2. 
Apparently the beneficial impact of the shadowing cannot 
be observed in this case, when the BS can be chosen only 
among two closest BS. A general remark is of the following 
order: strong shadowing requires larger geographical domain 
in which the serving BS is searched, as the optimal one may 
be located far from the mobile. 



C. Analysis of the path-loss factor 

We begin with an important remark regarding the scaling 
of E[?(0)] with respect to the density of the BS. 

Remark 4.5: Unlike the mean interference factor E[/(0)] 
(cf. Remark 4.2), the mean path-loss factor E[^(0)] depends 
on the intensity A of BS. By the homothetic invariance of our 
hexagonal and Poisson models, it is easy to see in the case 
of the distance-loss function (2) that this dependence has the 



following form E[Z(0)] = \'^/^ f E[/(0)] | a=i) . Consequently, 
in particular, the path-loss factor becomes preponderant in the 
case of sparse networks (small A) and negligible for dense 
networks (lai-ge A). We will see in Section V-B that E[Z(0)] 
can be evaluated explicitly in the case of the infinite Poisson 
network with an arbitrary distribution of the shadowing. 
Figures 5(a) and 5(b) show the mean path-loss factor E[Z(0)] 
evaluated for the intensity of BS A = 1.155BS/km^ (equivalent 
to A = 1km). The main observations are presented in the next 
section. 

D. Conclusions on numerical results 

For the hexagonal network we have observed the following 
facts regarding our two QoS "pre-metrics". 

• The mean path-loss factor increases to infinity in the stan- 
dard deviation of the shadowing, increases in the path- 
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Fig. 5. Mean path-loss factor expressed in dB, in hexagonal and Poisson network on the torus Tjv with log-normal shadowing with log-SD v and path- 
loss exponent /3. Note that E[i(0)] decreases with the size of the network. The analytical expression for the infinite (on R^) Poisson model is given in 
Proposition 5.6. 
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Fig. 6. Graphical explanation of a possible shape of the dependence of the 
blocking probability on the variance of the shadowing. 



loss exponent, increases in the cell radius, but (slightly) 
decreases in the number of base stations. 
• The mean interference factor is not monotone in the 
standard deviation of the shadowing: first increases and 
then decreases to 0. It decreases in the path-loss exponent, 
is invariant with respect to the cell radius and increases 
in the number of base stations. 

Knowing that the blocking probability depends on some com- 
bination of the path-loss and interference factors of users (cf. 
formulas (7) and (6)), and having observed that the mean 
values of these two factors have opposite monotonicities in 
the path-loss exponent, it is not surprising that the blocking 
probability is not monotone in the path-loss exponent; cf. the 
first observation of Remark 4. 1 . A similar argument explains 
a possible non-monotonicity of the blocking probability in the 
standard deviation of the shadowing; cf. the second observa- 
tion of Remark 4. 1 and the scheme on Figure 6. 

For the Poisson network we have observed the same tenden- 
cies of QoS "pre-metrics" as for hexagonal network mentioned 
above, except that the mean interference factor is monotone 
decreasing in the shadowing. Moreover, for large standard 
deviation of the shadowing, the "pre-metrics" of the Poisson 
network are very close to those of the hexagonal network. 



In the next section we will prove also that for the infinite 
Poisson network the distributions of our QoS "pre-metrics" do 
not depend on the shadowing and admit explicit formulas for 
their means. 

V. Mathematical results 

In this section we will state and prove some mathematical 
results regarding E[/(0)] and E[/(0)], which support and 
extend the numerical findings of Section IV. 

A. Toroidal models 

We begin by a simple observation regrading the log-normal 
distribution of the shadowing S with mean 1. Recall, it can 
be represented as S* = er'^ /2+aN wjjgj-g ]\[ jg (jjg standard 
Gaussian random variable. Thus, for any fixed e > we have 

Pr{5 > e} = Pr{A^ > ct/2 + (log e)/(j} "^ , 

which shows that the random variable S converges in prob- 
ability to when a (and hence v — crlO/ log 10) tends to 
infinity (and this even if £[5*] = 1!). From this, we have 
that the path-loss between any location y and any BS X, 
Lx{y) = L(\X — y\)/Sx{y), converges in probability to 
infinity when the variance of the shadowing increases. Conse- 
quently, for any finite network $ of base stations, the path-loss 
factor l{y) = min;^g$ Lx{y) converges in probability and in 
expectation to infinity. This explains the asymptotics of E[Z(0)] 
for large v observed on Figures 5(a) and 5(b). 

The somewhat surprising observation on Figures 4(a) 
and 4(b) regarding the beneficial impact of the strong log- 
SD V of the shadowing on the mean interference factor can 
be also confirmed mathematically. 

Proposition 5.1: Assume an arbitrary, fixed, finite pattern 
{Xi, X2, . . . , Xn\ ofBS locations. Consider any deterministic 
distance-loss function < L{r) < 00 and (independent) log- 
normal shadowing S'xi(') with the log-SD v. Then for any 
location y we have lim,j_j.oo /(?/) =0 in probability. 



Proof: It is enough to show liniu^oo Prj /(y) > e} = 
for any e satisfying < e < 1. Denote by Gi = 
Sxi{y)/L{\Xi — y\) the path-gain from Xi to y. Consider 
ordered vector (G(i), . . . , G(„)) of these path gains, where 
min^ Gi = G(i) < . . . < G(„) = maxj G^. Note that f{y) = 
l/^Cn) EiLi G(j) - 1 < {n ~ l)G(„_i)/G(„). In order to 
prove our claim it is enough to show that Pr{ G(„_i)/G(„) > 
e } — >■ when w — > oo. To this regard denote L{\Xi — y\) ~ li, 
and recall from the definition of our path-loss model that 
we can represent Gi{y) = e^\ where {Ni]i=i....^n are 
independent Gaussian random variables, with mean E[iVi] — 
— \ogli — cr^/2 and the same SD g — u log 10/10. Since 
Gi is monotone increasing in Ni we have G(i) — e^(*'. 



where min.^ Ni = iV(i) < 



< N. 



(n) 



maxi Ni. Moreover, 



A := { G(„_i)/G(„) > e } = { iV(„) - Ar(„_i) < M }, where 
M = -lege. Denote by Ay = {0 < Ni - Nj < M}. 
Note that A C IJ^ .^ „ ^ / Aij and the result follows 
from the fact that for any i y^ j Pr{Aij} — > when 
w — > oo. Indeed, for i ^ j, Ni — Nj = N is Gaussian 
random variable with mean log{lj/li) and variance a^ and thus 
Pt{A,j} = Pr{0 < TV < M} -> for any given finite M 
when a^ = v^ log 10/100 — > oo. This completes the proof. 

■ 

Corollary 5.2: Assume Poisson or hexagonal network on 
the torus Ttv, with log-normal shadowing having log-SD v. 
Then the mean interference factor /(O) converges in distribu- 
tion and in expectation to when w — > oo. 

Proof: For any e > 0, by Proposition 5.1 and Lebesgue 
dominated convergence theorem we have Pr{ /(O, $) > e } = 
E[Pr{/(0, 1") > e|$}] -)■ 0, when w ^- oo. This proves 
that /(O) converges in distribution to 0. Convergence of 
E[/(0)] to follows again from the Lebesgue dominated 
convergence theorem by the observation /(i/, $) < ^(Ttv) — 1 
and E[$(TAr)] < oo. ■ 

Remark 5.3: Recall that the log-normal distribution of the 
shadowing is heavy-tailed. The result of Proposition 5.1 and 
Corollary 5.2 can be heuristically explained and conjectured 
for other heavy tailed shadowing distributions by the so called 
single big jump principle, cf e.g. [^ , Section 3.1]. It says that 
the only significant way in which a large value of the sum of 
independent heavy-tailed variables can be attained is through 
a big value of single term of the sum ("big jump"). In other 
words, the maximum and the sum of independent heavy-tailed 
random variables have the same asymptotic of the distribution 
function for large values. Note also, that we observe the "single 
big jump principle" in a different scenario: we study the ratio 
of the interference (sum of the log-normal path-gains minus 
the largest path-gain) to the serving-BS path-gain (the largest 
one) asymptotically for large variance. 



B. Infinite models 

In this section we will consider infinite hexagonal and 
Poisson models. We will show first that serving BS X^, and 
hence the path-loss and interference factors, are well defined. 
Then we will argue that values of these factors in the infinite 
models can be seen as limits of respective toroidal models 
on Tm when n ^ oo. Finally we will prove a (surprising?) 



invariance of E[^(0)] and E[/(0)] in the infinite Poisson model 
with respect to the distribution of the shadowing. In this case 
the values E[^(0)] and E[/(0)] can be evaluated explicitly. 

Proposition 5.4: Consider infinite Poisson $ = ^p or 
hexagonal $ — $// model of BS, with shadowing whose 
marginal distribution has finite moment of order 2/ j3 (^). Then 
there exist X^ € <1> satisfying (3). Moreover, the path-loss 
factor and the interference factor calculated with respect to 
the restriction of $ to Tn, i-e., Z(0, $'"'") and /(O, $'"'"), 
converge almost surely and in expectation to Z(0,$) and 
/(O, <&), respectively. 

Proof: To prove the first statement it is enough to 
show that the expected number of BS Xi such that 
SxM/L{\X^\) > M is finite for any M < oo. In the 
case of the Poisson p.p. this will be shown in the proof of 
Proposition 5.6 below. Here we consider only hexagonal case 
$ = ^H- Denote by G{x) = Pr{ S" > a; }. We have 

n#{X^ e $H : 5x.(o)/i(|X,|) > M}] 
= e[ ;^ l(^Sx,(o)>MLi\X, 



< 



Xie<s>f 

e[ J2 g[ml{\x,\ 

oo 

J2QnG({nAK/2f/M) < oo, 



where 1() denotes the indicator function and the 
last inequality follows from the assumption E[S'^/^] = 
2//3 J„°° s'^/^~^G{s) ds < oo. This completes the proof of the 
first statement. 

In order to prove the second statement, note that for any 
realization the network $, for N large enough Xq e T^. 
Consequently, Z(0,$'^") is eventually constant in N while 
/(O, <&^") eventually increases in N (the serving BS is not 
changing any more and only interference is added). The 
convergence of the expectation of the path-loss factor fol- 
lows from the monotone convergence theorem, noting that 
l{0, cE>^") is decreasing in N. The convergence of the expec- 
tation of the interference factor follows form the dominated 
convergence theorem knowing that /(O, $) < /'(O, $), where 
/'(O, $) is the interference factor calculated under assumption 
that the handover policy selects the geographically closest BS 
as the serving one. By the independence of the shadowing 
fields given the locations of BS and the assumption that the 
mean shadowing is equal to 1 



E[/'(0,l>)]=E 



E 






L{\X'„*\) 



L X 



- 1. 



(8) 



where Xq* is a point of $ closest to the origin 0. By our 
assumption on the mean path-gain E[l/S'] < oo. The second 
expectation (8) is equal to the mean interference factor in the 
infinite model with constant shadowing S = I, and it is known 
to be finite in the infinite hexagonal and Poisson model; cf. 
respectively Remark 5.5 and Proposition 5.6 below. ■ 

^i.e., E[5^/'^] < oo. Note that 2//3 < 1 and thus the above assumption 
follows from our default assumption E[S] = 1 < oo. 
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Remark 5.5: It was shown in [8] that in the case of 5 = 1 
and the deterministic distance-loss function (2) E[?(0)] and 
E[/(0)] in the hexagonal model can be approximated by the 
following expressions 

^r./„ ^ M 0-9365 

E[/(0,$ff)] 



/3-2 



E[Z(0,$h)] 



K^ 



(7rA)/5/2(l + /3/2)- 

To the best of our knowledge, analytical expressions (approx- 
imations) in the case of the infinite hexagonal network with 
random shadowing are not known. We consider now infinite 
Poisson model. 

Proposition 5.6: Assume infinite Poisson network, deter- 
ministic distance-loss function (2) and a general distribution 
of the shadowing S satisfying £[5*^'^] < oo. Then the 
distribution of the interference factor /(O) = /(O, $) does 
not depend on the distribution S and the distribution of l(Q) 
depends on S only through the product XE[S'^'^]. Moreover 

E[/(0)] = ^ 



E[;(o)] 



/3-2' 



(ttAE [527/3] )/3/2 ' 

where T{a) ~ L t"^^e^*df. In particular, for the log-normal 
shadowing 

^,^,,n^l i^^r(l + /3/2)exp[(l-2//3)a2/2] 

Remark 5.7: The above result says that in the infinite 
Poisson notwork the existence of shadowing has no impact 
on the mean interference factor The impact of the shadowing 
on the mean path-loss factor in this model consists in a 
"fictitious" scaling of the intensity of the BS by the factor 
(E[S'2'^])^'2 < I xhe respective expressions in the case 
of 5 = 1 has been found for the first time (to the best 
of our knowledge) in [17]. Note however, that the above 
observation is valid only if the handover policy selects the 
BS with the smallest path-loss, as described in Section III-C. 
Indeed, assume that, despite non-constant shadowing, the 
handover policy selects the geographically closest BS as the 
serving one. Then, the mean interference factor E[/'(0)] can 
be expressed as in (8). Recall that the second expectation in 
this expression is equal to the mean interference factor in the 
same model without shadowing (i.e., S = 1). By the Jensen's 
inequality E[l/S'] > 1/E[S'] ~ 1 and consequently we 
observe the increase of the mean interference factor compared 
to the "shadowing-dependent" handover policy. In particular, 
for log-normal S with mean 1 and log-SD v we have 
E[l/S'] = e""' = e"'i°s'i"/i"o, which means that the log- 
normal shadowing in any geometric model of BS in which 
it is not taken into account in the handover policy increases 
the mean interference factor by u^ log 10/10 dB, where v is 
log-SD of the shadowing. 

Proof of Proposition 5.6: Note that the values of lifi) and 
/(O) are entirely defined by the collection of random variables 
{Lx{0) = L(|X|)/S'x(0) : X e $}. Given $ these random 
variables are independent. Thus by the displacement theorem 



for Poisson p.p. (cf. [23, Theorem 1.3.9]) {Lx(0)} = * 
constitutes a (non-homogeneous) Poisson p.p. on M+ = [0, oo) 
of intensity measure A' given by 

A'([0,.]) = E[vI/([0,.])] 



A / Vr{L{\z\)/S<s}dz 

/>oo 

27rA / rVv{L{r)/S < s) dr 
Jo 



27rA 



= 27rAE 



AS^/Z^TT, 



K^ 



E 



rE [1 (i(r)/5 < s)] dr 
r dr 



Sis 



Note that the latter expression is finite, which proves that 
the serving BS Xq is well defined (cf. proof of Proposi- 
tion 5.4). Note also that it depends on the shadowing only 
through its moment E[S'2/'^]. Moreover one obtains the same 
expression in the model without shadowing and the density 
of BS multiplied by E[S'2/^]. By the homothetic invariance 
of the Poisson model with the distance-loss function (2) the 
distribution of /(O) does not depend on the intensity of the 
BS. Thus the invariance of the distribution of /(O) on the 
distribution of the shadowing. In particular, we can conclude 
that E[/(0)] = 2/{j3 - 2) — the value obtained in the model 
without shadowing; see [17], cf. also [23, Example 4.5.1]. 
The formula for the mean path-loss factor follows from its 
dependence on the intensity of the base stations via the 
function X^^l"^ . This completes the proof. ■ 

VI. Concluding remarks 

We show that the QoS in path-loss-and-interference limited 
cellular networks is not always decreasing in the strength 
(variance) of the log-normal shadowing, provided the handover 
policy selects the BS with the smallest path-loss as the serving 
one. Under strong shadowing it principally suffers from the 
poor path-loss conditions with respect to the serving BS. 
For moderate shadowing however, when the QoS is not yet 
compromised by the path-loss conditions, it may profit from 
the reduction of the interference. This is because increasing 
variance of the log-normal shadowing tends to "separate" the 
strongest (serving BS) signal from all other signals — the 
phenomenon observed for heavy-tailed distributions and called 
"single big jump principle". This mathematical result seems 
also to be in line with a recent real-network observation [ ] 
that mobiles in indoor communications (typically subject to 
strong shadowing) report fewer BSs for potential handover 
The results presented in this paper regard the network-average 
of the QoS metrics. More study is needed, to analyze the 
impact of the shadowing on the distribution of these metrics 
in the network. This requires appropriate models of the spatial 
correlation of the shadowing. 
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